clear all; close all; 
k = sym(['960939379918958884971672962127852754715004339660129306651505519271702802395266424689642842174350' ... 
'718121267153782770623355993237280874144307891325963941337723487857735749823926629715517173716995' ... 
'165232890538221612403238855866184013235585136048828693337902491454229288667081096184496091705183' ... 
'454067827731551705405381627380967602565625016981482083418783163849115590225610003652351370343874' ... 
'461848378737238198224849863465033159410054974700593138339226497249461751545728366702369745461014' ... 
'655997933798537483143786841806593422227898388722980000748404719']); 
[x,y] = meshgrid(0:1:106, 0:1:16); 
tupperformula = rem(floor(floor((y+k)/17) .* 2.^(-17*x - rem((y+k), 17))), 2); 
tupperformula = double(tupperformula); 
image(fliplr((1-tupperformula)*255));
colormap gray
axis equal 
title('Tuppers Self Referential Formula'); 
set(gca,'YTick', []); 
ylabel('1.63 * 10^{49} \leq y \leq 1.62 * 10^{49} + 17'); 
